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We analytically study the effect of a uniform shear flow on the persistence properties of coarsening 
systems. The study is carried out within the anisotropic Ohta-Jasnow-Kawasaki (OJK) approxima- 
tion for a system with nonconserved scalar order parameter. We find that the persistence exponent 
9 has a non-trivial value: 9 = 0.5034 ... in space dimension d — 3, and 9 = 0.2406 . . . for d = 2, the 
latter being exactly twice the value found for the unsheared system in d — 1. We also find that the 
autocorrelation exponent A is affected by shear in d = 3 but not in d — 2. 
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I. INTRODUCTION 

The phenomenon of persistence in nonequilibrium sys- 
tems has attracted considerable interest in recent years 
EL both theoretically and experimentally 

H, 13 ■ The persistence probability, Pq (t) , of a fluctuating, 
spatially homogeneous nonequilibrium field is the prob- 
ability that the field X(t) at a given space point has not 
changed sign up to time t. This probability typically de- 
cays as a power law, Po(t) ~ t~ e at late times, where the 
persistence exponent 9 has in general a nontrivial value. 
Persistence has been studied in a considerable number 
of systems such as simple diffusion from random initial 
conditions, phase-ordering kinetics, fluctuating interfaces 
and reaction-diffusion processes 0. 

Experiments to determine persistence exponents have 
been carried out in the context of breath figures [f| , liq- 
uid crystals @, soap froths 0, diffusion of Xe gas in 
one dimension [8j and fluctuating monatomic steps on a 
metal/semiconductor adsorption system Si-Al surface 0. 
Many of these cases are examples of coarsening phenom- 
ena, where a characteristic length scale increases with 
time as the system relaxes towards an equililibrium that 
it attains only after infinite time in the thermodynamic 
limit. The experimental results are generally in good 
quantitative agreement with (exact or approximate) the- 
oretical predictions. 

A classic example of a coarsening phenomenon is the 
dynamics of phase ordering, where a system is quenched 
from a disordered high-temperature phase into an or- 
dered low-temperature phase. In the simplest case of a 
two-phase system, domains of the two equilibrium phases 
form and grow with time. The characteristic length scale 
at a given time is the typical scale of the domain structure 
that has formed at that time. The coarsening dynamics 
is usually characterised by a form of dynamical scaling, 
in which the system looks statistically similar at different 
times apart from an overall change of scale ^(J • 

Recently there has been interest in the effect of shear 
in a variety of systems such as macromolecules, binary 
fluids and self-assembled fluids Shear introduces 

anisotropy into the spatial structure. For systems under- 
going phase ordering in the presence of shear, the domain 
growth becomes anisotropic and this results in different 



growth exponents for the structure along and perpen- 
dicular to the flow. At present it is not clear whether 
shear leads to a stationary steady state, or whether do- 
main growth proceeds indefinitely at asymptotically large 
times |l2J. Shear may also induce phase transitions: for 
example, shear-induced shift of the phase transition tem- 
perature in the microphase separation of diblock copoly- 
mers has been observed Il3l. 



In this paper we analytically study the effect of an 
imposed uniform shear flow on persistence for the sim- 
plest case of a nonconserved scalar order parameter. We 
exploit a version of the Ohta-Jasnow-Kawasaki (OJK) 
approximation in phase-ordering kinetics ,10], modified 
to account for the anisotropy induced by the shear [l4| . 
Persistence is defined here as the probability that a point 
comoving with the flow has remained in the same phase 
up to time t. We employ an approach called the in- 
dependent interval approximation (IIA) which has been 
successfully used to obtain rather accurate values for per- 
sistence exponents in unsheared systems p] • This proce- 
dure assumes that the intervals between zeros of the pro- 
cess X(t) are statistically independent when measured 
in the mapped time variable T = Int. We find that the 
persistence exponent 9 is nontrivial and dimensionality 
dependent. For d = 3 we find 9 ~ 0.5034, compared to 
~ 0.2358 in the unsheared caseQ, while 9 ~ 0.2406 for 
d = 2 compared to 9 ~ 0.1862 without shear |3J. Re- 
markably, the value of 9 in d = 2 is exactly twice the 
value obtained for the unsheared system in d = 1 Q 
using similar methods. There is a technical subtlety in 
d = 2 which requires a careful definition of the persis- 
tence probability. In both d = 2 and d = 3 the shear 
increases the persistence exponent. 



The paper is organised as follows. In the next section 
the OJK theory is introduced and the autocorrelation 
function, which is a necessary input to the IIA calcula- 
tion, is obtained for d = 3 and d = 2. Section III contains 
a brief outline of the IIA, the results of which for the 
sheared problem are presented in section IV. Concluding 
remarks are given in section IV. 
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II. THE OJK THEORY 



In fc-space, Eq. Q can be written as 



We consider a nonconserved scalar order parame- 
ter <f>(x, i) evolving via the time-dependent Ginzburg- 
Landau equation 



d(j)(x, t) 
dt 



= V 2 <P(x,t)-V (</»), 



(1) 



where V((f) is a symmetric double- well potential. The as- 
sumption that the thickness £ of the interface separating 
the domains is much smaller than the size of the domains 
allows one to write an equation of motion for the inter- 
face, called the Allen-Cahn equation [l^. The velocity v 
of the interface is proportional to the local curvature and 
given by 



v{x, t) = —X7.n(x,t), 



(2) 



where n(x, t) is the unit vector normal to the interface, 
defined in the direction of increasing order parameter. 
The normal vector can be written in general as 



n(x, t) 



Vm(f, t) 
\7m(x,t) | 



(3) 



where m(a?, t) is the smooth field that has the same sign 
as the order parameter <f> and vanishes at the interfaces 
(where the order parameter vanishes). It is easier to work 
with an equation of motion for m(x,t) than for (f>(x 7 t), 
an idea that is exploited in the OJK theory 

By considering a frame locally comoving with the inter- 
face, with a space-uniform shear in the y-direction and 
flow in the x-direction (i.e. the fluid velocity profile is 
given by u = ^yex), where 7 is the constant shear rate 
and e x is the unit vector in the flow direction, the OJK 
equation for the field m(x,t) becomes [l4j 



dm(x, t) 
dt 

V 2 m(x,t) ~ 



iy- 



dm(x, t) 
dx 



lJab\P)- 



a,b=l 



dx a dx b 



where 



D a b{t) = (n a n b ) 



(4) 



(5) 



and (...) denotes average over initial conditions (or, equiv- 
alently, over space). The correct equation for m involves 
an unaveraged D ab , but the equation is then nonlinear 
and intractable. The essence of the OJK approximation 
is the replacement of the product n a rib by its average. For 
an isotropic system this gives, by symmetry, D a b = S a b/d, 
and the equation for m reduces to the diffusion equation. 
For the anisotropic sheared system, however, D a b{t) has 
to be determined self-consistently [Tj|- From Eq. JSJ, it 
follows that 



£a,«(*) = i. 



(6) 



dm(k,t) dm(k,t) 



Of 



dk,. 



- fc ° 2 + D ab (t)k a k b m(k, t). (7) 



a,fc=l 



The case d 



We now consider the above equation in dimension d = 
3 and solve it via the following change of variables , 

(k x ,k y ,k z ,t) -> (q x ,q y - jk x T,q z ,r), (8) 

with the introduction of an equivalent field 

fj,(q,r) =m(k,t). (9) 

The left hand side of Q now becomes dfi/dr and as 
a result equation (0 can be integrated directly to give 
(after transforming back to original variables) 



m(k, i) = m(k x ,k y + jk x t,k z ,0) 



x exp 



1 3 

-J k aM ab {t)h 



■ib=l 



(10) 



with non-vanishing matrix elements 

M n (t) = R 11 (t) + 2 1 tR 12 (t)+ 1 2 t 2 R 22 , 

M 12 (t) = Ri 2 {t)+itR 22 {t), 

M 22 {t) = R 22 (t), 

M 33 (t) = R 33 (t), 

where 



(11) 



Rn(t) 


= 4 


1 dt'{[l-D 11 (t')]+2 1 t'D 12 (t') 
j 2 t' 2 [l-D 22 (t')}}, 




+ 


Rl2(t) 


= 4 


[ df{-D 12 (f)-7*'[l--D22(f)]} 
Jo 




= 4 


[ df{l- D 22 (t')}, 


R 3 3(t) 


= 4 


[ dt'{l-D 33 (t')}. 



(12) 



Due to the symmetry of the original OJK equation 
the terms Ri 3 , R 23 , M\ 3 and M 23 all vanish. The as- 
sumption that the initial condition, m(fc,0), has a Gaus- 
sian distribution, appropriate to a quench from the high- 
temperature phase, is used throughout the paper. 

In order to use the IIA to investigate the persistence 
properties of the coarsening system, it is first neces- 
sary 0, Q to compute the autocorrelation function of 
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the rescaled field X(t) = m(x, t)/{[m(x, t)} 2 ) 1 / 2 , which is 
constructed to have unit variance, using the initial cor- 
relator 



The quantity ( 
give 

(lm(x,t)} 2 ) 



m(x, 0)m(x', 0)) = AS d {x ~ x'). (13) 
m(x, t)} 2 ) 1 / 2 can easily be evaluated to 



1/2 



(2tt)3/ 2 ^Det M(t) 



1/2 



(14) 



Turning now to the two-time correlator of X(t), we 
recall that we want to calculate this correlator not at a 
fixed point in space, but at a point that is advected with 
the shear flow. Due to the shear, the field at the space- 
time point (x + jyti,y, z,t\) at time t\ will be at the 
space-time point (x + 72/i 2 , y, z, t 2 ) at time t 2 . The auto- 
correlation function a(ti,t 2 ) = (X(ti)X(t 2 )) is therefore 
given by 



(2tt) 3 , 

^— -f- ^Det M (ti)Det M(t 
A 2 



1/2 



The next step is to evaluate the term 
/c-space implies 



(m(x + jyti,y, z, ti)m(x + jyt 2 ,y, z, t 2 )) . (15) 
in the above equation. We note that average over initial conditions in 



(m(k x ,k y + jk x h, k z , 0)m(k' x , k' y + jk' x t 2 , k' z ,0)) = (2n) d A5(k x + k' x )5{k y + 7M1 + K + jk' x t 2 )S(k z + k' z ). (16) 



Using Eqs. (jl()|l andl^EJ) we can evaluate the term 



i(l)m(2)) = A^exp 

k 

A 



1 3 
-- 2] k a B ab k b 



a, 6=1 
1 



(2tt)3/2 ^Dct B(t u t 2 y 



(17) 



where 
Bu(t) 

B 12 (t) 
B 22 (t) 
B 33 (t) 
Bi 3 (t) 



M u (ti) + M n (t 2 ) + j 2 M 22 (t 2 ) x (i a - ti) 2 

-2 7 {t 2 -h)M 12 (t 2 )]/2, 

[M 12 {h) + M 12 (t 2 ) - 7 (t 2 - t x ) M 22 (t 2 )}/2, 

[M 22 {ti) + M 22 (t 2 )}/2, 

[M 33 {h) + M 33 (t 2 )]/2, 

B 23 {t) = 0. (18) 



The notation (1) and (2) in (|17l) denotes space-time 
points (x + jyti, y, z, ti) and (x + ~fyt 2 ,y, Z,t 2 ) respec- 
tively. The problem is now reduced to evaluating the 
determinants of the matrices M(t) and B(tx,t 2 ), as the 
autocorrelation function o(ti, t 2 ) can now be written as 



a(t\,t 2 ) 



[DetM(ti)DetM(i 2 )] 
vTJet5(ti,i 2 ) 



1/4 



(19) 



The terms M a t,(t) cannot be computed explicitly for gen- 
eral t; only in the scaling limit (i.e. t — > 00) can one make 
progress. In this limit it can be shown that (to leading 
order for large t) |l4| 



Mn(i) 



15 



7 2 t 3 , 



Mu(t) 



It 1 



M 22 (t) = -t, 
16 

M 33 (t) = —t. 

5 



(20) 



Using equations l)2()|l. the determinants of M(i) and 
B(t\,t 2 ) can now be evaluated leading to 



Dct M(t) 
DetB(*i,t 2 ) 



375 ' ' 
8 7 2 (*2 + 



\5 r 



125 
1 - 



2U 2 



2t 2 H x 

(t 2 +hr 



(t 2 +h) 2 

The autocorrelation function follows: 

1 w 5 /\t u t 2 ) 



(21) 



a(h,t 2 ) 



2 [l-fW'fr.t*)] 1 / 2 ' 
1 sech 5/2 (T/2) 



[l-fscch 4 (T/2)] 



1/2 • 



(22) 



where W(h,t 2 ) = Mit 2 /(t 2 + ti) 2 , T = T 1 -T 2 and the 
final form follows after introducing the new time variable 
Ti = hit;. After this change of the time variable, the 
autocorrelation function depends only on the time differ- 
ence Ti — T 2 . Since the process X(T) is also gaussian, 
the process X(T) is a gaussian stationary process. This 
will be the case whenever the autocorrelation function of 
X depends on t\ and t 2 only through the ratio t\/t 2 , i.e. 
when it exhibits a scaling form. 
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B. The case d = 2 

For d = 2, we follow the same analysis as for d = 3 but 
with the change of variables 



*) (q x ,q y - -fk x T,T). 



(23) 



The solution of Eq. J7J in the large-* limit is given by the 
d — 2 analogue of Eq. lfTU|) : 



m{k,t) = m(k x ,k y + jk x t,0) 



x exp 



I 2 

-- 22 k a M ab (t)h 



a, 6=1 



(24) 



The matrix elements M a b(t) can be evaluated for large f 
using asymptotic analysis along the lines outlined in ref. 
H, with the result 



3 7 * 2 



Mu(t) = 47* 2 v / ln7*- 



Mi 2 (t) = At^/h^t- 



2* 



VEt* 



Afaa(t) = -VM* 5 

7 



(25) 



where we have retained just the leading subdominant 
terms, of relative order l/ln(7*). 

The subleading terms in Mn(i) and Mi^if) are nec- 
essary as there are cancellations to leading terms in the 
determinant of M(i), which is given by DetM(*) = 4* 2 . 
Using Eq. (|2"Hl the following averages can be calculated: 



([m(x,t)ff 2 = 



1 



(m(l)m(2)) 



A 



A_ 

2tt ^/DetM(t) 
1 



1 1/2 



(2tt) ^Det B(h,t 2 )' 



(26) 



where the matrix elements B a b are given by the expres- 
sions in Eq. (|18|) but with the corresponding M a t,(t) given 
their by their d = 2 equivalents in Eq. (|25|l . The autocor- 
relation function a(t\,t 2 ) for <i = 2 can now be evaluated 
using the set of equations (|2"6*)l to give 



a(*i, * 2 ) 



1 1/2 



4ti* 



1'2 



*?( 1+ \/S) +i 2 (l + V /§|) 



(27) 



Note that a(ti, t%) given by that Eq. (|27|l does not have 
a scaling form, i.e. it is not simply a function of *i/*2j due 
to the logarithms. However it does a scaling regime. In 
the limit t\ — > oo, * 2 — > with ti/t 2 fixed but arbi- 
trary, the ratio of logarithms can be replaced by unity 
and a(ti,t 2 ) depends only on t\/t 2 in this regime. In 



terms of the new time variable T 
tains 



ln(*2/*i), one ob- 



a(ti,t 2 ) = v / sech(T), 



(28) 



where T — T\ — T 2 , i.e. the process X(T) becomes sta- 
tionary in the defined scaling limit. We will use Eq. I|28|) 
rather than Eq. I|27|l to extract 9 for d = 2, but one 
must note the special limit taken to derive (|28|l where the 
persistence probbaility P{t\, t 2 ) is the probability that a 
point moving with the flow has stayed in the same phase 
between times t\ and t 2 . 



III. THE INDEPENDENT INTERVAL 
APPROXIMATION 

The above analysis in both d = 3 and d = 2 shows that 
X(t) is stationary in the new time variable T (with the 
caveat noted above for d = 2). We note that the expected 
form for the probability, Po(t), of X(t) having no zeros 
between t\ and t 2 , namely Po ~ {t\/t 2 ) 9 for t 2 ^> *i, 
becomes exponential decay Po ~ e - 6 '( T 2- T i) i n the new 
time variable. This reduces the problem of calculating 
the persistence exponent to the calculation of the decay 
rates |l7| . 

The order parameter field in the OJK theory is given 
by (f> — sgn(X). The autocorrelation function 



A(T) = (mm) = (sgnX(O)sgnX(T)), 



(29) 



for the field at a space point moving with the flow, is 
given by 



A(T) = - sin" 1 a(T), 

7T 



(30) 



which follows from the fact that d is a Gaussian field 
We will determine the persistence probability Po (t) from 
A(T). 

We briefly discuss the IIA |l| and use it to obtain ap- 
proximate values for the exponent 6 following the de- 
velopment in ref. 0. In the scaling limit, the interfaces 
occupy a very small volume fraction and as a result 4>(T) 
takes values ±1 almost everywhere. The correlator A(T) 
can be written as 



A(T) = J2(-l) n Pn(T), 



(31) 



n=0 



where P n (T) is the probability that the interval T con- 
tains n zeros of 4>{T). For n > 1, P n (T) is approximated 
by, assuming that the intervals between zeros of X are 
independent, 

P n (t) - (T)- 1 f dTt [ dT 2 ■ ■ ■ / dT n 

^Q(T X )P(T 2 - Ti) • • ■ P(T n - T„_i)Q(T - T n ), 

(32) 
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where (T) is the mean interval size, P(T) is the distribu- 
tion of intervals between successive zeros and Q(T) is the 
probability that an interval of size T to the right or left 
of a zero contains no further zeros, i.e. P(T) = —Q'(T) 
where the prime indicates a derivative. The II A has been 
made in Eq. (|32|) by writing the joint distribution of zero- 
crossing intervals as the product of the distribution of 
single intervals. The Laplace transform of Eq. <\','<2\) leads 
to P(s) = [2 - F{s)]/F(s) where 

n*) = 1 + ^r 5 ! 1 - sA (s)} (33) 

and A(s) is the Laplace transform of A(T). 

It is straightforward to show that the mean interval size 
is (T) = -2/A'(0). The expectation that P (T) ~ e~ 6T 
for large T implies a simple pole in P{s) at s = —9. The 
persistence exponent 9 is therefore given by the first zero 
on the negative axis of the function 

s T 2s f°° 
F(s) = 1 - — [l - - 1 dT exp(- S T) 

X sin -1 o(T)] . (34) 

For further analysis is is useful to first extract the 
asymptotic behavior of the autocorrelation function a(T) 
of the field X(T). From Eqs. JZZ) and (J22J) we find, for 

T —>■ oo, 

We now turn to the results. 



IV. RESULTS 

The term A'(0) can easily be evaluated to give A'(0) = 
-y/17/2/n in d = 3 and — v/2/tt in d = 2. From $Q 
F(Q) = 1, and from iLiol) F(s) diverges to — oo for s — > 
—5/4 and —1/2 in d = 3 and 2 respectively. Therefore, 
the zero of F(s) lies in the interval (-5/4,0) and (-1/2,0) 
for d = 3 and 2 respectively. Solving 13411 numerically 
for this zero, we get the IIA values for the persistence 
exponent as 9 = 0.5034 . . . for d = 3 and 9 = 0.2406 . . . 
in d = 2. In the absence of shear the IIA gives [3| 6 = 
0.2358 . . . in d = 3 and 0.1862 . . . in d = 2, which agree 
quite well with simulations 19]. 

A very interesting feature of the d = 2 result for 9 
is that it is exactly twice the value of the exponent ob- 
tained within the same approximation (i.e. using OJK 
theory and the IIA) for the unsheared problem in one 
space dimension: 0f^ 2 = 20*^- That this must be so 
is easily seen directly from the form (|28|l for a(ti,t 2 ) for 
the sheared problem in = 2. The equivalent result 
for the unsheared system in general space dimension is 
a(t 1 ,t 2 ) = sech d/2 (T/2) 3]. For d = 1 this is identi- 
cal to Eq. H28fl apart from an overall factor 2 in the 



(logarithmic) timescale T. It follows that the relation 
®sh 2 = 2#un S h d° es not require the IIA but only that the 
underlying field m (or, equivalently, X) be Gaussian, i.e. 
it requires use of the OJK theory but not the IIA. It is 
interesting to speculate that it might even hold beyond 
the OJK approximation, in which case one might imag- 
ine that there is a very simple explanation for it. As yet, 
however, we have been unable to find one. 

The autocorrelation function A(t\,t 2 ) is also interest- 
ing. In the limit t 2 3> t\ that defines the autocorrelation 
exponent A [l^, via A(ti,t 2 ) ~ (tx/t 2 ) x , the quantity 
a(ti, t 2 ) is small and Eq. (|30|l can be linearised in a(ti, t 2 ) 
to give, from Eq. fJSJ with T = ln(i 2 /fi), 

A(t ^~{fc/S'\ "d = 2. <36) 

These results give A = 5/4 for the sheared system in 
d = 3, compared to A = 3/4 in the unsheared system 
|10|. whereas for d = 2 the autocorrelation exponent 
takes the same value, A = 1/2, in both cases. We should 
repeat the caveat that, for d = 2, the simple power-law 
form H36|) requires the limit t± — > 00, t 2 — > 00 with t 2 /ti 
fixed but large. If t 2 — > 00 for fixed t\, Ea. (|27|l gives 
a(h,t 2 ) - (ii/t 2 ) 1/2 [ln(7t 2 )/ln(7ti)] 1/4 , which does not 
have a simple scaling form (it is not simply a function of 

h/t 2 ). 



V. CONCLUSION 

We have studied the effect of shear flow on the persis- 
tence exponent 9, for a system with nonconserved scalar 
order parameter, using an approximate analytical ap- 
proach based on the OJK theory and exploiting the "in- 
dependent interval approximation". The persistence is 
defined in a frame locally moving with the flow. 

The exponent 9 is nontrivial and is increased by the 
presence of shear. This implies that the shear accelerates 
the change of sign of the fluctuating field. In dimension 
d — 2 we find the intriguing result that 9 has a value equal 
to twice that of the unsheared system in d = 1 , within the 
OJK theory. The autocorrelation exponent A increases in 
the presence of shear for d = 3 but is unchanged by the 
shear in d = 2. 

For nonconserved dynamics in the absence of shear, 
experiments on liquid crystals have been performed to 
measure both 9 H and A [2(j . There is also a recent ex- 
periment on the measurement of a two-time correlation 
function in order-disorder phase transition in CU3AU |2l| . 
Liquid crystal experiments are a possible candidate for 
testing our predictions in a model with nonconserved or- 
der parameter, and numerical simulations may also prove 
useful. On the analytical front, the method of the cor- 
relator expansion jj] might be used to obtain a more 
accurate result for 9 in d — 3 than can be obtained using 
the IIA. 
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